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E fiber elastic modulus, MPa 

e”’ matrix elastic modulus, MPa 

Egff effective modulus of the matrix in the loading direction, MPa 

Ej^ unloading elastic modulus at N number of cycles 

Eq initial elastic modulus of the first cycle (modulus of undamaged 

laminate), MPa 
E^ secant modulus 

Esds laminate modulus assuming damaged matrix material, MPa 

R stress ratio, S j /S 

' * min' max 

Smav maximum laminate stress, MPa 

max 

minimum laminate stress, MPa 
S ^1 laminate stress in the 0° fiber direction, MPa 

Sgg laminate stress in the 90° fiber direction, MPa 

Y cyclic-hardened yield stress, corresponds to one-half the matrix 

fatigue limit for R = 0, MPa 
Ae laminate strain range 

AeJ’oii^P compressive strain range of the matrix material in the loading 
direction 

AS laminate stress range, MPa 

AS$^ stress range that causes no fatigue damage, MPa 
0 ^ axial stress in the fiber in the loading direction 

axial stress in the matrix in the loading direction 
stress in the matrix material in the loading direction when the 
matrix yields 
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ABSTRACT 

Previous research has shown that boron/aluininum can develop significant 
Internal matrix cracking when fatigued. These matrix cracks can result 
in a 40 percent secant modulus loss in some laminates even when fatigued below 
the fatigue limit. The present study shows that the same amount of fatigue 
damage will develop during stress or strain controlled tests. Stacking 
sequence has little influence on secant modulus loss. The secant modulus loss 
in unidirectional composites is small, whereas the losses are substantial in 
laminates containing off-axis plies. This paper presents a simple analysis 
that predicts unnotched laminate secant modulus loss due to fatigue. The 
analysis is based upon the elastic modulus and Poisson's ratio of the fiber 
and matrix, fiber volume fraction, fiber orientations, and the cyclic-hardened 
yield stress of the matrix material. Excellent agreement was achieved 
between model predictions and experimental results. With this model, designers 
can project the material stiffness loss for design load or strain levels and 
assess the feasibility of its use in stiffness critica'' parts. 

INTRODUCTION 

Metal matrix composites (MMC), in spite of their relatively high cost, 
have several inherent properties that make them attractive for structural 
applications; MMC have high stiffness- to-weight ratios, high strength-to- 
wel ght ratios, and better transverse strength, better operative temperature range, 
and environmental resistance than do competitive epoxy-resin-matrix composites. 
Many components are of continuous fiber-reinforced MMC, such as boron/aluminum 
(B/Al). 
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Previous research [1-3] has shown that boron/aluminum can develop significant 
internal matrix cracking even when cycled below the fatigue limit. This 
results in laminate modulus loss. In quasi-isotropic laminates, matrix cracks 
reduce stiffness as much as 40 percent. Because most MMC structural components 
are expected to be stiffness critical, even a small drop in component stiffness 
may render the part useless or cause failure of the structure. 

This paper builds upon the data base developed in References [1-3,9] by 
adding rejiUlts for three additional laminates. These additional tests are 
conducted under both strain and load control. References [1-3,9] dealt with 
fatigue mechanisms while this paper produces information that is more design 
applicable. The damage model analysis presented in Reference [1] is modified 
and expanded to predict the laminate secant modulus after approximately 
500,000 fatigue cycles (when the laminate is in a saturation damage state). 

The analysis needs only the elastic modulus and Poisson's ratio of the fiber 
and matrix, fiber volume fraction, fiber orientations, and the cyclic-hardened 

w 

yield stress of the matrix material. With this analysis, designers can 
project the material stiffness loss for design load or strain levels and 
assess the feasibility of its use in stiffness critical parts. 

The secant modulus is predicted for eight different laminates and compared 
to experimental test results, The tests were conducted at numerous stress 
and strain levels below the fatigue limit. Test and predictive results are 
discussed. 


EXPERIMENTAL PROCEDURE 

The material tested was unnotched boron- aluminum composites with a 
6061 aluminum matrix and 0.14 nni diameter boron fibers. Table 1 presents 
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material properties for the boron and aluminum constituents and Table 2 shows 
the eight laminates that were tested, Notice that five laminates ([Ojg, 

[0/90]25» [90/0]25» C0/i4{V90/0/i45/'§Qr]g and [0/i45/903g) are from previous 
references and three laminates (C0/±45]g, and [±45/023g) are new. 

The laminate specimen widths are shown in Table 2, All specimens were annealed 
such that the 6061 aluminum matrix was in the annealed condition before 
testing. All specimens were fatigue loaded at 10 cycles per second except 
when the stress-strain response of the material was recorded on an x-y plotter. 
The stress-strain data were taken under quasi -static condition periodically 
during test life. The strain was measured with a 25.4 mm gage- length 
extensometer. All specimens tested for this study were cycled at constant 
amplitude stress or strain levels below their fatigue limit. The stress 
ratio, R, was constant for each test presented. Different specimens were 
tested at R ranging from 0,0 to 0.5, 

To obtain the data generated in Reference [3], fatigue tests were con- 
ducted at a constant amplitude stress range until failure or to two million 
cycles. The secant modulus of approximately the 500,000th cycle was obtained 
from recorded stress-strain data. Therefore, each tested specimen provided 
one data point fcr the 500,000th cycle response. The tests conducted in the 
present study (laminates [Og/i^Sjg, [MB/Ogjg, and [0/±45]g) were somewhat 
different. To generate more data per specimen, tests were conducted at a 
constant stress or strain level for 500,000 cycles (time enough foir a 
saturation damage state to develop) and then the stress-strain response was 
recorded. The stress or strain range was then increased to a new desired 
level, and another 500,000 cycles applied. The resulting stress-strain 
response was recorded. This process was repeated up to as many as five 
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different stress or strain levels per specimen. Fatigue damage at each level 
depends only on the applied stress range and therefore is not influenced by 
the prior cycling at lower stress ranges [1], 

BACKGROUND 

If fatigue damage is to be avoided in general, and low cycle fatigue 
failures in particular, the cyclic loading must produce only elastic strains 
in the constituents. Even so, local plastic straining can be permitted in 
the composite during the first few load cycles, provided that the composite 
"shakes down" during these few cycles. The shakedown state is reached if 
the matrix cyclically hardens to a cyclic yield stress Y such that, 
subsequently, only elastic deformation occurs under load cycles. The 
shakedown limit for the composite containing 0° fibers is considerably 
below the composite's fatigue limit. Previous tests have shown that the 
matrix fatigue limit coincides with the stable cyclic yield stress for 
annealed aluminum [3,4,5] and [5]. The value of Y is 70.38 MPa for 
annealed 6061 aluminum [2,3]. 

The shakedown stress range for a unidirectionally loaded laminate can 
be found by using laminate theory to determine the yield surfaces for 
individual plies in the laminate. Figure 1 shows an example of a [0/±45/90]g 
layup under biaxial inplane stresses S-ji and $ 22 * Each ply has its own 
elliptical yield surface, constructed analytically from the ply matrix 
stresses and the von Mises yield condition. The overall yield surface of the 
laminate is an internal envelope of the yield surfaces of individual plies. 

The shakedown stress range, ASg^, is the width. 6f the overall yield surface 
in the direction of uniaxial loading applied in this experimental program. 


The value of can be easily calculated with computer program AGLPLY [6], 

More details of this procedure can be found in References [3,7]. 

When a specimen was cycled above its shakedown range* matrix cracks 
were observed [3]. Figure 2 shows a micrograph of cracks in a 45“ ply 
matrix. These matrix cracks reduce the effective tensile modulus of the 
matrix. The cracks tend to open and close under remotely applied cyclic 
loads. This matrix cracking and the subsequent crack opening and closing 
results in a bilinear response as will be explained later in the paper, but 
may be observed from experimental stress-strain nisponses as shown in 
Figure 3 for the 500,000th cycle. The amount of damage (matrix cracking) 
can be inferred from the changes in elastic unloading modulus.j as 
described in References [1-3]. Stress-strain data were taken at intervals 
during the fatigue cycling to record the drop in laminate modulus as a 
function of the number of cycles. The damage was expressed in terms of its 
effect on normalized by E^. An example of the fatigue damage accumula- 
tion as a function of number of applied cycles and stress level is presented 
in Figure 4 for a [0/A45/90/0/i45/'5S^g laminate. Most damage occurred 
in the first 500,000 cycles. Notice that each specimen appears to reach a 
stabilized value of E^j/Eq. herein referred to as a "saturation damage state" 
(SOS). After the saturation damage state is reached, the laminate will 
neither accumulate more damage nor fail under the present loading condition. 

Returning to Figure 3, the cyclic stress-strain curve for the fourth cycle 
results from elastic and plastic defornntion and has a secant modulus, E^, 
of 10.82 N 10^ MPa, For the same specimen, the 500,000th cyclic stress-strain 
curve has a very different shape with an associated secant modulus of 
8.88 X 10^ MPa. The shape change of the stress-strain curve and the drop in 
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secant modulus (almost 22 percent) Is attributed primarily to matrix cracking. 

In contrast, hardening of the matrix material [2,3] usually caused the secant 
modulus to increase after some initial cycling. Notice that the fourth cycle 
in Figure 3 has a secant modulus much less than the elastic modulus (Table 2) 
just due to plasticity. If the laminate was cycled at or below the shakedown 
range, the matrix would harden such that the secant modulus would be approximately 
equal to the elastic modulus. To confirm that the decrease in elastic 
unloading modulus and secant modulus observed during fatigue cycle was 
caused by matrix cracking and not fiber breakage, the aluminum matrix was 
gradually etched in a 30 percent HCA solution in distilled water. Fiber 
failure was detected only in specimenb tested at stresses near the fatigue 
limit. However, substantial laminate modulus changes were detected for 
specimens stressed well below this level. Specimens that sustained modulus 
loss had long matrix cracks that grew parallel to the fibers in the off-axis 
layers of the laminate. Some cracking perpendicular to the loading direction 
has been observed in the matrix of the 0“ plies [8,9]. (Notice in Figure 1 
that the 45° and 90° plies yield at a lower laminate stress than do the 
0° pliesi therefore, the matrix in the off-axis plies would be expected to 
undergo more plastic deformation which, in turn, would lead to more fatigue 
cracking than in the 0° plies. Indeed, the off-axis plies did have more 
cracks.) The individual cracks did not extend into adjacent plies that had 
different ply orientation. No delamination was discovered between the plies, 
as is commonly reported for polymer matrix composites. Therefore, practically 
all of the observed elastic unloading modulus decrease was caused by cracks 
in the matrix, since such cracks were the only observed damage of consequence. 
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DAMAGE MODEL 

A version of this damage model was first presented in Reference [1,3], 
The present model is quite similar except that the bounds on the matrix 
stress range, Aa"’, are modified. Whereas the previous model was intended 
to calculate only laminate unloading elastic modulus, the present model is 
extended to calculate the secant modulus, E,, 

A simple analysis is developed to relate the decrease in laminate secant 
modulus caused by matrix damage. The model starts with the matrix cycling 
plastically, As cracks develop due to plastic cycling, the effective modulus 
is reduced for the portion of the matrix cycle that is in tension. The 
model presents simple equations to approximate the effective matrix modulus 
due to the cracking at an assumed cyclic strain range. The program AGLPLY 
is used to calculate the laminate response with the effective modulus of the 
fatigued matrix. Thus the bilinear response illustrated in Figure 3 is 
computed. The secant modulus is calculated from the bilinear response. 

Figure 5 illustrates this behavior in terms of the applied laminate 
stress and the corresponding axial stresses in the matrix and 0° fibers. The 
dashed lines in Figure 5 represent the initial loading response. Accordingly, 
the first load cycle causes the matrix and 0® fiber stresses to follow the 
dashed loops. The laminate has an ideally elastic-plastic matrix (for 
illustration of the model and simplicity of presentation) and is subjected 
to a constant cyclic stress range, AS. The dashed loops are for the same 
condition represented in Figure 2 in the fourth cycle, is assumed to be 

the axial stress in the matrix material in the loading direction at the 
shakedown stress limit, ASg^, (the matrix is yielded at this point by a 
combination of axial and shear stresses). Assuming the matrix yields at the 
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same value in tension and compression, equals half of the laminate's 
shakedown strain range, ASgj^/E^, times the inatri) tensile modulus, e”’. 



lET * ^ 


( 1 ) 


The in this equation is the shakedown stress range, E^ is the undamaged 

laminate’s elastic modulus in the loading direction, and e"' is the undamaged 
matrix’s elastic modulus, (og^ is approximately equal to Y for unidirec- 
tion composites. The axial matrix stresses at yielding was assumed to be 
Y and -Y in the previous model [1],) With subsequent cycling, the cyclic 
plasticity causes matrix cracks to initiate and grow, effectively decreasing 
the matrix tensile modulus until a saturation damage state is reached. The 
dashed loops in Figure 5 narrow to zero-width loops, shown as solid lines, 
which represent the saturation damage state. These solid lines correspond to 
the laminate cyclic stress-strain response illustrated in Figure 3 for the 
500,000th cycle. The saturation damage state develops when the matrix 
cracking causes the load to transfer to the 0° fibers, thus relieving the 
matrix from undergoing additional damaging plastic deformation. 

The drop in matrix modulus in the load direction due to fatigue damage 
will now be determined using Figure 6. The strain in the matrix and laminate 
is plotted versus the matrix stress, o”’, and laminate stress, S, respectively. 
The damage state has an associated cyclic strain range, At. If this cyclic 
strain range is assumed, an effective tensile modulus e”’^^ of the matrix 
material can be estimated. This assumes thatthe same SOS will be reached 
by either stress or strain control. Note that is the modulus in the 

loading (0“ fiber) direction. The compressive strain range of the matrix, 
vvas approximated by 
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( 2 ) 




m 

comp 


m 


“sh 


The effective tensile modulus of the matrix material can now be approximated 
by dividing Og^ by the cyclic strain minus the compressive portion. 



m 


°Sh 


Ac - Ae 


m 

comp 


(3) 


The is used as the matrix modulus in lamination theory (using 

the computer program AGLPuy [6]} to calculate the unloading elastic modulus 
of the composite in its saturation damage state, Egpg, (at approximately 
500,000 cycles). The shear modulus of the matrix is also reduced within 
AN6PLY based on and Poissor»*s ratio. All the fibers were assumed to 

be intact, the matrix damage was assumed to be characterized by the laminate's 
lower modulus, Although such a formulation implicitly assumes that 

the matrix modulus is reduced isotropically, the reduction really is orthotropic. 
However, interest is in the laminate modulus in the primary loading direction 
only, and the assumption should not introduce excessive error, 

Returning to Figure 6, we now know the modulus for each of the two 
linear segments, as well as the strain ranges. Therefore, the overall 
laminate stress range, AS, can be calculated as follows: 

AS - (A^comp^^o " ^^comp^^SDS* 
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Equatton (4) was rewritten using Equation (2), 


* ^SDS "S’ ^ 

* Eq ^ ^hh’ 

The values of ASg|^» Eq» and Egpg came from AQLPLY* Equation (4) applies 
for either stress or strain control cycling, By selecting a number of different 
strain range values, Ac, the corresponding laminate stress range, AS, 
can be calculated and plotted versus Ac. The laminate secant modulus Is 

Eg * AS/Ae. (6) 

RESULTS AND DISCUSSIONS 

The predicted cyclic stress-strain response after 500,000 cycles and 
Its associated secant modulus is presented In this section and compared 
with measured experimental results. The predictions are shown as solid 
lines (see Figure 7 as an example), For reference a dashed line representing 
the undamaged elastic modulus of the laminate is shown. The secant modulus 
scale can be read in two ways. First, by entering on the AS axis, crossing 
to the solid prediction line and down to the secant modulus scale; this would 
give the secant modulus of a laminate after 500,000 cycles at a given stress 
range. Secondly, one could simply drop from the cyclic strain scale directly 
to the secant modulus scale to assess the secant modulus after 500,000 cycles 
at a given strain range. Notice that the secant modulus scale is nonlinear. 

Also notice that the secant modulus scale ends on the left at the shakedown 
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limit; the secant modulus is equal to below the shakedown limit. The 
experimental data were generated at stress ratios, R, between 0.0 and 0.5. 
Since all the data show little scatter, this confirms that the damage 
developed in the matrix is a function of stress range, AS, and not of R 
(or mean stress) [3]. 

Figures 7 through 12 present the experimental and analytical correlation 
in order of elastic nwdulus, highest to lowest. 

Figure 9 and 10 presents data for the [ 05 /± 45 ]g and [0/±45]g laminates, 
respectively. These tests are significant because they were conducted under 
both stress and strain control. The experimental data indicate that the 
same damage state is reached whether the stress is held constant and the 
strain increases or the strain is held constant and the stress decreases. 

This material behavior allows one to assume a constant strain range to 
calculate fatigue darnage for strain or stress control tests in the 
presented analysis. 

Figure 8 and 11 include data points representing the Initial cyclic 
response of the laminate (e.g., the fourth cycle in Figure 2). These data 
illustrate the secant modulus loss due to matrlH yielding. The Initial 
cyclic stress-strain responses are reasonably close to the predicted response 
after 500,000 cycles, however* as shown in Figure 3, the reason, or mechanism 
for the secant modulus loss is different. 

Figures 8 and 9 show that stacking sequence has very little effect on 
the secant modulus, Eg, in L‘0/90]2g - [90/0]2g, and [0g/±45]g - C±45/02]s 
laminates, respectively. Previous research [1] showed that stacking 
sequence may have an effect on the degradation of the elastic unloading 
modulus, E(^, in particular, near the shakedown limit. 
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In general, the predictions fit the experimental data very well. The 
data fell slightly above the predictions In some cases and slightly below 
Ifi others. Some of this scatter may be attributed to deviation of the fiber 
volume fraction from what was assumed. Also, the annealing treatment of the 
aluminum matrix may have varied from one laminate to another. It Is also 
acknowledged that the assumption of Isotropically decreasing the modulus 
of the matrix due to cracking may affect the predictive results for various 
laminates differently. In any case, the present model does a very good 
job of representing the extent of accumulated fatigue damage in the saturation 
damage state and predicting the obss ved material response. 

The experimental data fit the predictions quite well, even though the 
individual data points were generated at different stress ratios. This 
confirms observations [3] that the matrix damage Is a function of stress 
range and not mean stress. 

The resulting secant modulus after 500,000 cyles Is significantly below 
the elastic modulus for all of the tested laminates, except the [0]g 
laminate. If compared at a cyclic strain range of 0.004, the [0]g laminate 
retained approximately 95 percent of the original elastic modulus. The 
other laminates retained about 60 to 70 percent of their original moduli. 

These differences between the often calculated elastic modulus and the 
resulting secant modulus must be addressed by the designers of stiffness 
critical parts. Certainly the unidirectional laminate may still retain the 
desired stiffness, but laminates with off-axis plies must be scrutinized as 
to their design load levels and stiffness requirements. 
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SUMMARY 

A simple model has been presented that will accurately predict the 
secant modulus loss, under fatigue loading, of boron/a lumlnum laminated 
composites. The model requires only knowledge about the fiber and matrix 
nwduli and Poisson's ratio, fiber volume fraction, layup orientation, and 
the cyclic hardened yield strength for annealed materials (the fatigue 
limit of the matrix material). These properties, in conjunction with lamination 
theory and von Misses yield criterion are sufficient to predict modulus 
changes. 

Since the assumptions are not restricted to boron fibers and no empirical 
constants are used, the author feels that this model can be applied to other 
continuous fiber reinforced ntetal matrix composite systems, particularly 
those with the matrix in the near annealed condition. 

The experimental data indicated that the same degree of damage can be 
reached from either strain or stress control testing; this verifies a basic 
assumption in the model. The experimental data further indicated that laminate 
stacking sequence did not have a large influence on secant modulus drop. 

The presented results indicate that laminates with off-axis plies may result 
in secant moduli that are about 25 to 40 percent below the elastic moduli 
after fatigue cycling. 

Hopefully, this model for metal matrix composites will help designers 
more fully understand these materials and their limitations. It is further 
intended to help the designers avoid the pitfalls that are caused by generating 
S-N fatigue data alone, without regard to the behavior of the material below 
the fatigue limit where the data in this paper were generated. 
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TABLE 1.- COMPOSITE CONSTITUENTS MECHANICAL PROPERTIES 





TABLE 2.- DIMENSIONS AND PROPERTIES OF BORON/ALUMINUM LAMINATES 
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Figure 1.- Yield surfaces for plies of a B/Al laminate. 
The $22 direction coincides with 0° fiber 
direction. 
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Figure 5.- Decreasing tensile matrix response 
and Increasing 0“ fiber stress as the 
laminate attains a saturation domage state. 



Figure 6,- Composite laminate and matrix stress-strain 
response for a saturation damage state. 
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Figure 8.- Correlation of experimental and model predictions for 
CO/90123 and [90/0123 after 500^000 fatigue cycles. The 

cyclic response for the ^tth cycle is shown for some specimens. 
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Figure 9.- Correlation between experimental and model predictions for 

[02A^5] and [-^5/02]^ laminates after 500^000 fatigue cycles. 
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Figure 10,- Correlation between experimental and model predictions 
for [0/±45]o laminate after 500,000 fatigue cycles. 
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